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Abstract: We study the effects of the temperature and of a magnetic field in the setup of an 
intersection of D3/D7 branes, where a large number of D7 branes is smeared in the transverse 
directions to allow for a perturbative solution in a backreaction parameter. The magnetic 
field sources an anisotropy in the plasma, and we investigate its physical consequences for the 
thermodynamics and energy loss of particles probing the system. In particular we comment 
on the stress-energy tensor of the plasma, the propagation of sound in the directions parallel 
and orthogonal to the magnetic field, the drag force of a quark moving through the medium 
and jet quenching. 
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1 Introduction 

One of the amazing developments emerging from the research in string theory, is the idea of 
a gauge/gravity correspondence [1]. The remarkable feature of this correspondence is that it 
relates the strongly coupled regime of the gauge theory to the weakly coupled regime of the 
string theory and vice- versa. Consequently, it has become a powerful tool in studying strongly 
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interacting systems by using a conjectured dual weakly coupled string/gravitational theory. 
At present, holographic descriptions of non-perturbative phenomena include, among others 
applications to condensed matter physics, high energy physics and quark- gluon plasma. 

One of the most distinctive uses of the gauge/gravity correspondence has been the study 
of the physics of heavy ion collisions. Through collisions at Brookhaven and LHC a strongly 
coupled plasma of quarks and gluons was created which cannot be described by the standard 
perturbative techniques. Also other methods such as Lattice Gauge theory fail in com- 
puting transport coefficients of the plasma and the rapid thermalization rate of the quark 
gluon plasma observed. This is where the gauge/gravity duality enters in the field and pro- 
vides interesting new insights. For example, for large-iV gauge theories at strong coupling, 
gauge/gravity duality predicts that the ratio of the shear viscosity to entropy density is 1/4-7T 
(see [2]) in natural units and therefore very close to the measured value. The small value 
of the ratio of the shear viscosity to entropy density can be understood as an effect of the 
strong coupling of the system. Within the framework of Gauge/ Gravity duality we can also 
compute the thermalization rate Tth ~ 0.5fm of the plasma [3, 4], which is in agreement with 
the observed value indicating again the strong coupling nature of the plasma. The energy 
loss of the heavy/energetic partons in the plasma also acquires a gravity dual description (see 
[5])- 

Despite the remarkable insights into the quark-gluon plasma and QCD in general gained 
by studying gauge/gravity dualities, the application of the correspondence to real- world sys- 
tems such as QCD remains a challenge and has to be developed further. So far, we do not 
have a rigorous string dual of QCD at hand. However, under extreme external parameters 
(such as temperature and chemical potential) different gauge theories exhibit similar proper- 
ties. Therefore, it is natural to apply holographic techniques to study phenomena which are 
believed to be of universal nature. 

An important example in this class of phenomena is the effect of mass generation and 
spontaneous chiral symmetry breaking in the presence of an external magnetic field. The 
effect is known as magnetic catalysis and has been shown insensitive to the microscopic 
physics underlying the low energy effective theory. Using conventional field theory methods, 
the magnetic catalysis has been demonstrated in various (1+2) and (l-f-3)-dimensional field 
theories [6] , while the holographic study of the effect initiated in [7] 1 . Additional holographic 
studies of magnetic catalysis at finite temperature or chemical potential appear in [9]. 

Until recently all the holographic studies of the magnetic catalysis were in the probe ap- 
proximation, where the backreaction of the flavor branes on the supergravity background is 
neglected [10]. On the field theory side, this corresponds to an approximation in which the 
flavor degrees of freedom Nf are much smaller than the color ones N c . Unquenching the holo- 
graphic description means a large number of flavor branes that backreact on the geometry. 

1 For a comprehensive review we refer the reader to [8]. 
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Due to the technical difficulties that arise from a set of localized flavor branes, we distribute 
them along the compact directions [11]. This procedure is called smearing 2 and restores a 
significant part of the global symmetry of the geometry. 

A promising framework for the construction of such a geometry was started in [14], where 
the ten-dimensional supergravity solutions including the backreaction of a large number of 
D7-branes in AdSs x (with X§ any squashed Sasaki-Einstein manifold) was introduced. 
This was further developed in [15], where the black-hole solution dual to the non-conformal 
plasma of flavored M = 4 super symmetric Yang-Mills theory is presented 3 . The authors 
outline the smearing procedure, derive the corresponding equations of motion and present a 
perturbative solution for general massless non-supersymmetric flavor D7-brane embeddings. 

The first steps towards unquenching the holographic description of magnetic catalysis have 
been undertaken in [19] and [20]. More specifically in [19], a string dual to SU(iV c ) M = 4 
SYM coupled to Nf massless fundamental flavors in the presence of an external magnetic field 
is presented. For sufficiently strong magnetic field, the supergravity background is unstable, 
suggesting that the theory undergoes a phase transition to a stable phase with dynamically 
generated mass for the matter fields. In [20], the external magnetic field couples to Nf massive 
fundamental flavors and the background has a hollow cavity in the bulk of the geometry, 
where it is similar to the supergravity dual of a N = 1 non-commutative SYM. The radius 
of this cavity is related to the dynamically generated mass of the fundamental fields. After 
developing an appropriate renormalization scheme, the free energy and the condensate can be 
expanded in powers of the perturbative parameter. While at leading order, both agree with 
the previously obtained results in the probe approximation, at next to leading order the effect 
of magnetic catalysis is enhanced and the contribution to the condensate runs logarithmically 
with the finite cutoff Auv- 

An overview of the paper is as follows: In section 2 we continue the studies initiated in 
[19, 20] and present a string dual to the finite temperature SU(iV c ) M = 4 SYM coupled to 
Nf massless fundamental matter in the presence of an external magnetic field. The solution 
is analytic and perturbative in a parameter that counts the number of internal fundamental 
loops. Given the non illuminating expressions for the functions of the background we provide 
some numerical plots, and since we have a perturbative solution we supplement it with a 
hierarchy of scales. 

In section 3 we study the thermodynamics of the anisotropic black hole, which provides a 
non trivial check for the validity of the gravity solution. Since the solution is first order in 
the expansion parameter, our computations have some overlap with those of [21] and extend 
those of [15] in the presence of an external magnetic field. While in the absence of a magnetic 

2 For a detailed review on the smearing see the review [12], while for other solutions employing this technique 
that appeared after the review see [13]. 

3 All the hydrodynamic transport coefficients of the model were analyzed in [16], while the addition of a 
finite baryon density was presented in [17]. For a review on unquenching the Quark Gluon Plasma see [18]. 
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field the breaking of conformal invariance happens at second order in the expansion parameter 
[15, 17, 22], in its presence conformal invariance breaks at first order. 

In section 4 we holographically calculate the stress energy tensor of the boundary field 
theory. The presence of the magnetic field sources an anisotropy in the medium, which is 
realized through a difference between the pressure transverse to the magnetic field and the 
pressure along the direction of the magnetic field. We present thermodynamic arguments 
supporting the holographic computation. 

In section 5 we calculate the energy loss of the partons as they propagate through the 
anisotropic plasma. The jet quenching parameter depends on the relative orientation between 
the anisotropic direction, the direction of motion of the parton and the direction along which 
the momentum broadening is measured. We consider a parton moving parallel to the magnetic 
field with the momentum broadening taking place in the transverse plane. The presence of the 
magnetic field enhances or reduces the jet quenching parameter of a theory without magnetic 
field, depending on the conditions we use to make the comparison. The drag force experienced 
by an infinitely massive quark propagating at a general angle through the plasma is calculated 
using an appropriate set up to compensate the Lorentz force on the probe quark. In this way 
we obtain an expression reflecting the anisotropy of the plasma due to the external magnetic 
field. 



2 Constructing the black hole 

The present section is devoted to the construction of a supergravity background describing 
an anisotropic black hole. The field theory duals are realized on the intersection between a 
set of N c color D3-branes and a set of Nf, homogeneously smeared, flavor D7-branes, with 
an additional coupling between the fundamental fields and an external magnetic field. 

2.1 Setup 

The smearing of the flavor D7-branes allows for an ansatz where all the functions of the 
background depend just on the radial coordinate. Having this in mind and inspired by 
[15, 19, 20], we adopt the following ansatz for the metric 



b 2 T dt 2 + b (dx\ + dxf) + dx 



+ h,2 



b 2 b\ S° F z do 1 + S' z ds z cp2 + F' z (dr + A CP 2) 



(2.1) 



where the CP 2 metric is given by 



ds 



cp2 = -^d>x 2 + - cos 2 ^{d6 2 + sin 2 9dip 2 ) + - cos 2 ^ sin 2 ^{dtp + cos 6d(p) 2 



A CP 2 = — cos 2 ^ {dip + cos 6 dip) 



(2.2) 
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The range of the angles is < (%, 9) < n, < (ip, r) < 2tt, < tp < Air. The ansatz for the 
NSNS and the RR field strengths is given by 

B 2 = Hdx 1 A dx 2 , 2 = JdtA dx 3 , 
F 5 = Q c (1 + *)e(5 5 ) , F 1 = Qf{dr + A CP 2) , F 3 = d0 2 + 5 2 A iq , (2.3) 

where £(S§) is the volume element of the internal space 4 and Q c , Qf are related to the number 
of different colors and flavors in the following way 

Q c Vol(X SE ) Ar AQ f Vol(X SE ) 



Nr. 



N 



(2^ 5s a'2 ' ^ Vol(X 3 )g s " {2A) 

In our case Xse = S 5 and the X$ = S* 3 , a 3-sphere with volume 2tt 2 . The fact that the 
flavors are massless is encoded in the independence of F\ on a, see [14, 23] . All the functions 
that appear in the ansatz, h, br, b, S, F, J and H , depend on the radial variable a only. 
In the convention we follow, S and F have dimensions of length, b, bx, h, J and H are 
dimensionless and a has a dimension of length -4 . The function b in the ansatz for the metric 
reflects the breaking of the SO(l,3) Lorentz symmetry down to 50(1,1) x 50(2). The 
blackening function &t allows for the existence of solutions with a black brane, whose horizon 
sits at a position ah such that 6r(c/i) = 0, and which allows to study the field theory at finite 
temperature. 

Solving the lOd equation of motion for F$, we need to impose the following relation 



J' = Qc- 



5 bl 



h 



(H-H 



(2.5) 



where Hq is an integration constant. In the next subsection we will keep the function J and 
will see how this relation appears from an effective one-dimensional Lagrangian. 

2.2 Effective actions and equations of motion 

The action for the Type IIB supergravity plus the contribution from the Nf D7-branes in 
the Einstein frame is 

5 = S H B + Sfi , (2.6) 
where the relevant terms of the Sub action are 



Sub 



1 



2k io 



d w xyf^g 



R 



-d M $d M <S> - -e 2rS> F 2 - --e^F 2 

2 3! ' ' 



2 
1 1 



m 



2k\ q 



2 3!" " ( 3 ) 

and the action for the flavor D7-branes takes the usual DBI+WZ form 



1 



ll F 2 
25T( 5 ) 



C 4 AH 3 A F 3 



(2.7) 



Sfi = -TrJ2 



d 8 xe*V-det(G + e-*/ 2 J-) 



Or + C 6 A T 



(2- 



4 With / e{S 5 ) = Vol(S 5 ) 
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with T = I?2 + 27ra'.F. In those expressions B2 denotes a non-constant NSNS potential which 
will model the magnetic field, F the worldvolume gauge field and the hat refers to the pull- 
back of the quantities, along the worldvolume directions of the D7-brane. The gravitational 
constant and D7~brane tension, in terms of string parameters, are 

1 - Tl - 1 

2^ ~7s~ WM^ 1 ' (2,9) 

We plug our ansatze, (2.1) and (2.3), into (2.6) and integrate out all the directions except the 
radial one, since the dependence is trivial. After an integration by parts to get rid of second 
derivatives we obtain the following expression 

S eff = ^^ [ Zeffdo (2.10) 

where Vi^ is the volume of the Minkowski space and the one-dimensional effective lagrangian 
C e ff is given appendix A. Since the function J enters in the effective action only via its radial 
derivative, there is a first integration given by a conserved quantity. We fix this constant of 
motion in the following way 

^ =-Q c H => J' = e ~*Q°% {H - H ) . (2.11) 

which is precisely (2.5). The next step is to use (2.11) to eliminate J' in favor of H in (2.10), 
after performing the following Legendre transformation 

L-eff - L eff J 



SJ' 



(2.12) 

J'=J'(H) 

and then calculate the Euler-Lagrange equations from the transformed action (2.12). The 
equations of motion are given in appendix A. 

Setting Qf = in the transformed action, the Euler-Lagrange equations imply that a 
solution with H 7^ is given by (black) AdS^ x X$e with $ = and H = Hq constants. 
We will use this solution later on as a starting point to obtain a black brane solution with 
backreacted flavor in the presence of a non trivial H. 

It is worth noting that by demanding dj'£ e ff = —Q C H exactly, with H the value of the 
magnetic field in the unflavored limit, we are enforcing the field J to vanish when Nf — > 0. 
As such, J reflects magnetic effects by providing a field connected holographically to the 
magnetization of the system, as we will see. 

The equation for the blackening factor (A. 3) decouples from the rest and can be solved 
analytically 

&2 = e -K° , (2.13) 

where is a non-extremality parameter coming from the integration constants. The position 
of the horizon is at a — > 00, whereas the boundary would be at a = (there is an additional 
integration constant corresponding to a shift in a, which we set to zero). 
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Reduced five-dimensional action For the calculation of the stress-energy tensor in sec- 
tion 4, we find convenient to write as well a truncated five-dimensional action, obtained after 
integrating out the compact Sasaki-Einstein manifold in (2.1). Denoting the effective metric 
as g^u, the action is 



&5d — TT~2 I d J X- v /— g [Ckin + £pot] + — 2 / d 5 x£ t0 p , 

ZKg J ZK 5 J 

where the kinetic, potential and topological terms are given by 

Ckin = R[9] - y^/^/ - 20d^w - Id^d"* - ^e*-™tF iafp F>°"> 



(2.14) 



(2.15) 



12 



_$_20 

-e 3 



/ TT TTfiVp 

11 aup 11 > 



c 



pot 



, 4e f f+2w ^10w 



6) 



Sl e f f-8 W +2$ 



Q 



l e *-V-*» Biu/B V> 



Qc 
2 



C 



top 



4Q /e 

Q 



|+2/+2u. 



(2.16) 
(2.17) 



with the convention £ tx v zr = l for the completely antisymmetric symbol. To make contact 
with the ansatz presented in section 2.1 we identify k 2 = k 2 /Vse and 



f 



log 



w 



log 



H 3 = dB 2 , F 3 = dC 2 



(2.18) 



g^dx^dx 1. 



— 10 

e 3 



b 2 T dt 2 + b (dxl + dx 2 2 ) + dxl + e~ wf h~h 2 b 2 T da 2 



B2 = \^B nv dx^ A dx u = H(a)dx A dy , C 2 = ^C^dx? A dx u = J(a)dt A dz . 

This effective 5d action is not enough to study perturbations, though, since the truncation of 
fields that cancel in the specific background we are considering is not a consistent one [24]. 

2.3 Perturbative solution 

The system (2.11) and (A.4)-(A.9) allows for a systematic expansion of all the functions in 
power series of Qj, as defined in equation (2.4). In fact physically it is more relevant to 
expand in the parameter, e* 

e* = Qje^* , (2.19) 

which takes into account the running of the effective 't Hooft coupling (through the dilaton 
factor e *). We consider the following first order expansion in e* 

i? 4 



b = 1 + e*6i + 0{el) , h - 

S = r (l + e^ + Oie 2 )) , F 
$ = + + 0{e 2 ) , H 



{l + e*h 1 + 0(el)) , 
r (l + ejk F 1+ C>( e 2)) , 
F (l + e^ + Oiel)) . 



(2.20) 
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where R 4 = Q c /4- We define the new radial coordinate r, in such a way that the zeroth order 
expansion in of h becomes R 4 /r 4 

e -K° = l _ -A . (2.21) 

The extremal limit corresponds to sending the horizon radius to zero. It is also convenient 
to define the following parameter 

r 4 m = e-*-H*# , (2.22) 

The result is a coupled system of second order differential equations which can be decoupled 
by the transformations 

Ai = S 1 - F x , Ti = 4fi + 16 Si + 5/h , Ai = h x - h . (2.23) 

This allows us to write 

¥? + - *&L 9l = ^ll±^lk (2 24) 

1 r(r*-r^) Vl - rf 1 (r* - r*)^T^ ' ( } 



~„ r 4 + 3r* 16r 2 - _ 4r 4 



where 



= {6, A, T, A, $} , C{6,a,t,a,$} = {0, 8, 8, 3, 0} , 
^{6,A,r,A,$} = {0, 0, -16, -1,4} , B {6 ,a,t,a,*} = {"4, 2, -6, -1, 2} . (2.26) 

The solution to these equations of motion is described in appendix B. Let us comment here on 
the boundary conditions we impose. In our solution there are four scales. We have already 
introduced the first three: is the radius of the horizon and we impose the fields to be 
regular there; r m is associated to the magnetic field, and r* denotes the point at which we 
pierce the dilaton, this is, «3>(r) = + <fi(r) with <fi(r*) = 0. With this scale we defined e* 
and its interpretation is given in terms of the scale at which the gauge coupling is defined, 
since [15] 

1 a Nt , 

* = to 9 ' N ° e w (2 - 27) 

The fourth scale (which we will define as r s ) is the scale at which we paste the thermal 
solution presented in the appendix B to the T = (supersymmetric) one [15], i.e. we impose 
the following conditions 

h(r s ) = H x {r 8 ) = Ai(r s ) = , Ti(r a ) = ^ , A 1 (r s ) = ^ . (2.28) 

Notice that &T( r s) ^ 1, which is the supersymmetric solution. This is not a problem in 
Euclidean signature, since it can be solved by fixing the periodicity of the Euclidean time in 
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the solution without temperature, and we will use this in the following to compare the energy 
and free energy of both solutions. In Lorentzian signature it introduces an error of order 
(r/j/r s ) 4 , which is small provided <C r s . 

Prom now on we set r s — > oo, which corresponds to push the Landau pole to infinity, or 
more physically, to focus only in the IR properties of the theory. A UV completion of the 
system is not known even in the supersymmetric case. The following results can be understood 
as the leading terms in an r/j/r s expansion. At the same time, we will take = r^, therefore 
describing the value of the dilaton relative to its value at the horizon, which implies that the 't 
Hooft coupling A/j = 4:7rg s N c e® h is evaluated at the energy scale marked by the temperature. 

For completeness, let us mention that from (2.11) and (2.25) we have at first order in 



d r J 



1 



-d r 



rdrHx 



(2.29) 



Qualitative behavior of the solution Given the gargantuan form of the solution to our 
system at first order in e^, which can be found in appendix B, we give in this section a 
description of the different functions presented above. In this section some numeric work is 
presented, but in the rest of the paper we will restrict to analytic results. 

The function b± is easy to describe by focusing in its radial derivative, given by 



2rt 



r ( r 4 _ r 4) 



log 



+ \A 4 + < 



(2.30) 



which for r > and real non-vanishing r m is always negative (it is exactly vanishing if 
r m = 0), and asymptotes — > at large radius. As the boundary condition used in the 
integration is b\(r s ) = 0, we conclude that this function is a monotonically decreasing function 
of r for finite r m (exactly zero if r m = 0) with the maximum value at the horizon. 
Similarly, we can analyze the radial gradient of the dilaton correction 



ir 2 r 2 + / (r 4 + r 4 i) ( r 4 +r 4j 



r r 2 



(2.31) 



which is strictly positive for r > and real r m . In this case the boundary condition used 
to integrate the solution is </>i(r*) = 0, where r* will be identified eventually with the hori- 
zon position as the IR scale of our effective solution. At large radius the dilaton diverges 
logarithmically, signaling the presence of a Landau pole, as discussed in [15]. 

For the other functions present in our solution -namely Ai, Ti, Ai and H\— the gradient 
does not take a simple form that is worth writing, so we provide plots of the functions for 
several values of the parameters. For example, for h\ one has that, numerically, the radial 
gradient is strictly non-negative (zero if r m = 0), and Ai(r s ) = from the boundary condition, 
in a similar situation to the function b\ but with different sign for the gradient. In figure 1 
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we plot this quantity as a function of r/r^ for several values of r m /r^ = 0, 2, 5, 10 and observe 
that it has non- negative gradient, and approaches Ax — > as r — > r s (with r s — > oo in the 
figure) . 



1.5 



.2D 2:5 3:0" 3X".--"-4"-0 



Figure 1. Ai as a function of r/rh for several values of r m /rh = (blue straight line), 2 (orange 
dashed line), 5 (brown dotted line) and 10 (purple dotdashed line). To produce this plot the limit 
r s — > oo has been taken analytically first. 



As opposed to the previously presented cases, function Ti presents some structure. To 
start with, the boundary condition at r = r s changes and is given by T(r s ) = 2/9. However, 
when one works in the r s — > oo limit this boundary condition is modified to Ti(oo) = 1/2, 
which is the value of the function when r m = 0. For small values of the magnetic field scale 
(weighted by the horizon radius), r m /r/j < 1.23144, the value of Ti at the horizon is less than 
1/2, and after that specific value of the magnetic scale it is always larger than 1/2. We plot 
this behavior in figure 2. Given the analyticity of the function there is a minimum which, 
numerically, we determined to be at r m ~ 0.961 122r/ l . We have not found any characteristic 



0.49 




Figure 2. Ti(r/,) as a function of r m /rh- We observe a minimum at r m = 0.961122rh with value 
Ti = 0.485816 and the curve crosses Ti = \ again at r m = 1.23144r^. To produce this plot the limit 
r s —> oo has been taken analytically first. 



signature of the presence of this minimum of Ti(r^) in the plasma. 



In figure 3 we plot several examples of Ti as a function of the radial variable in three 
graphs, classified according to the value of the function at the horizon. All the curves present 
a minimum (on the horizon when r m < 0.961122r/j and on the bulk otherwise) and asymptote 
the r m = value (Ti = 1/2) at large radius. 




L5 2.0 2.5 3.0 3.5 40 f ' 
(c) r m > 1.23144r h 



Figure 3. Ti as a function of r/ru for several values of r m /rh represented by a blue straight line, an 
orange dashed line, a brown dotted line, and a purple dotdashed line, with values given respectively 
by (a) 0, 0.32, 0.64, 0.961122, (b) 0.961122, 1.05, 1.14, 1.23144 and (c) 1.23144, 2, 3.5, 5. To produce 
this plot the limit r s — > oo has been taken analytically first. 

We have not given an analytic expression for Ai because we couldn't find an easy way 
to write it, since it involves integrals of Legendre functions. However, from integrating the 
equation numerically we find that its behavior is very similar to that of b\ or Ai (with reversed 
sign), and we simply report here figure 4. We are not going to need to evaluate Ai anywhere 
in this work. The reason is that this is the mode describing the squashing in the compact 
Sasaki-Einstein manifold, but from the point of view of the 5-dimensional system it is just a 
scalar that does not enter explicitly in the 5-dimensional metric (its influence would be felt 
just via the equations of motion, but recall we have defined Ai precisely to decouple them). 
In this paper we will focus on the thermodynamics, stress-energy tensor and energy-loss of 
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probes in the system, which do not depend explicitly in the matter content of our theory, just 
in the 5-dimensional metric. 

t 




1.5 2.0 2.5 3.0 3.5 4.0 



Figure 4. Ai as a function of r/77, for several values of r TO /r/j = (blue straight line), 2 (orange 
dashed line), 5 (brown dotted line) and 10 (purple dotdashed line). To produce this plot r s has been 
taken to r s = 80rh in the numerics. We have checked that this value for r s gives indistinguishable 
results from those in figures 1 and 3. 



Finally, we present the flavor correction to the NSNS 2-form H±. As usual we take 
the r s — > oo limit analytically and we find, as it was the case for Ti, that the boundary 
condition is not Hi(oo) = but Hi(oo) = —1/4. One might be puzzled by the fact that 
the correction is not vanishing independently of the value of r m , in concrete when r m = 0, 
but recall that this correction is modulated by the flavorless value of the NSNS field strength 
H ~ + ChHi + 0(e/i) 2 ), therefore at vanishing magnetic field we have H = 0. As the 
value of the magnetic field is increased the correction gets smaller and smaller as can be seen 
in figure 5. 



Hi 



-0.05 
-0.10 
-0.15 
-0.20 
-0.25 



\T.5 2.-0—- 3.0 



3.0 3.5 4.0 



Figure 5. H% as a function of r/r^ for several values of r m /r/, = (blue straight line), 2 (orange 
dashed line), 5 (brown dotted line) and 10 (purple dotdashed line). To produce this plot the limit 
r s — > oo has been taken analytically first. 
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2.4 Hierarchy of scales and regime of validity of the supergravity solution 

The perturbative solution, that we present in full detail in the appendix B, needs to be 
supplemented with a hierarchy of energy scales (in terms of radial scales). Our analysis follows 
closely similar sections of [15, 19, 20], whose arguments we repeat here for completeness. 

As usual, for the Taylor expansions in (2.20) to be valid in the region < r < r s we need 
to separate the scale r s from the scale introduced by the solution to 4>\{r), which diverges 
logarithmically at large values of the radius r > r s , r s <C r/ l e 1 / eh . The requirement that we 
discard corrections in rh/r s implies that our perturbative corrections are much larger than 
the terms we discard, therefore e/j 3> r^/r s . Joining these two conditions we have 

e -iAfc < r ± < eh (2.32) 
r s 

which for large r s 3> (implying that the UV completion this theory needs is far from the 
IR, where we study the physical properties of the system) implies that 

< e h ~ \ h ^± < 1 . (2.33) 

The scale r m is associated with the magnetic field and can be arbitrarily close to r m = 0. 
For large values of r m (large magnetic fields/magnetization of the system, as we will see in 
the next section), requiring that our solution remains in the perturbative level sets up a top 
value. As can be seen from the plots given previously, the maximum value of the functions 
appearing in the solution is at the horizon, and from the asymptotic values given in appendix 
B it is easy to see that, at large values of r m , all the functions diverge at the horizon as r^/r^. 
Therefore we must impose 1 3> £h^m/ r h' which gives the condition 

\r m \ < ^ ■ (2.34) 
e h 

Similarly to [15], validity of the supergravity approximation requires to ignore closed 
string loops (N c 1) and a' corrections (A^ 3> 1), where is the effective 't Hooft coupling 
at the energy scale set by the temperature. In addition, validity of the smearing approximation 
suggests a dense distribution of flavor D7-branes. In summary we have 

{JV C ,JV>}»1, A h »l, eh = ^L«l. (2.35) 

Finally requiring that a corrections, which scale as X h , are sub- leading relative to flavor 
corrections, controlled by e^, requires 

\~h' 2 < e h . (2.36) 
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3 Thermodynamics 



In the previous section we presented in full detail the solution of an anisotropic black hole 
and now we will extract its thermodynamic properties. This will provide a non trivial validity 
check of the solution itself through the closure of the standard thermodynamical formulae. As 
in [15, 17], all quantities are obtained in power series of the perturbative expansion parameter 
and, therefore, the relevant thermodynamic relations are verified up to the relevant order. 



3.1 Smarr formula 

The temperature of the black hole is computed after imposing regularity of the Euclidean 
action. A simple computation using (B.16), (B.19) and (B.23) gives 5 



rh 



R 2 



l + \e h (3Ai - Ti - h) 



irR 2 



J r=r h 



1 + \ th ( 1 




(3.1) 



The entropy density is proportional to A%, the volume at the horizon of the eight dimensional 
part of the space orthogonal to the t, r plane (where t is the Euclidean time), divided by the 
infinite constant volume of the 3d space directions V3. Another simple computation using 
(B.16), (B.19) and (B.23) gives 



2n_As = Nhl 
n 2 V 3 27Ti?6 



(3Aj 



Ti 



N 2 c rj 
2irR 6 



r=r h 




N 2 7T 2 T 3 



1 + 



1 




(3.2) 



Note also combining (3.1) and (3.2) that 

sT 



2vr 



r 4 
r h 



(3.3) 



'10 



In principle this result is perturbative in eh and valid to order e\ in the present case, however, 
it is not difficult to show that the statement is true, independently of the expansion parameter. 

We define now the magnetic quantities. One natural identification for the magnetic field, 
B, is given by the value of the H field at the boundary, which from (2.22) is 6 



D 



r 2 R~ 2 



(3.4) 



Looking at (2.11), we see that J' and Hq are conjugate variables. The existence of the 
holographic duality implies that, if we associate Hq with the magnetic field then J has to 



5 In this section all quantities have corrections coming from C(e|) terms as well as ^t, where we are setting 
r a —> 00. 

6 Notice that we cancel a factor of e* fe by passing between the string and Einstein frames. 
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determine the magnetization density A4. This relation reads 



M = 



1 



SS, 



eff 



V 3 J 5H 



2k\ q 



J'dr 



N 2 
c AT 



27T 2 R A ~"reg , 



(3.5) 



and we will regularize the finite temperature result subtracting the zero temperature one. 
Using (2.29), to obtain the integral of J', we arrive to the following expression for the mag- 
netization 

N 2 C QfB 



M 



2tt 2 R a 2 



log 



r 2 h + \/ri + ri 



(3.6) 



The next step in the determination of the Smarr formula is the calculation of the internal 
energy. Starting from the ADM energy we have 



Eadm = \-V~9tt [ (fix^det g 8 {K T - K ) . 



(3.7) 



The eight-dimensional integral is taken over a constant time, constant radius hypersurface. 
The symbols Kt and Kq are the extrinsic curvatures of the eight-dimensional subspace within 
the nine-dimensional (constant time) space, at finite and zero temperature, respectively. Using 
the explicit solution in appendix A we have 



eadm 



Eadm _ 3iV c 2 r^ 



8vr 2 R 8 



l + ^log 



(3. 



Another way to write (3.8) is, at order 0{e 2 l ), 

3 



&ADM 



sT 



1 



BM 



(3.9) 



4 2 

which implies that we must identify the ADM mass with the magnetic enthalpy of the system, 
tt — £adm ■ The internal energy, U, is given in terms of the enthalpy by the following 
expression 

U = H + BM = ysT + \bM . (3.10) 



3.2 Thermodynamic potentials 

The relations that must be satisfied by the thermodynamic potentials are the following 

F = U-sT, Q = T - B M , (3.11) 

where J- is the Helmholtz free energy (in the ensemble where the magnetization is kept fixed) 
and Q is the Gibbs free energy (in the ensemble where the magnetic field is kept fixed), which 
is the interesting ensemble in our case. These thermodynamic potentials are related by a 
Legendre transformation 



T 



dT 
dM 



M . 



(3.12) 
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In a holographic set-up the thermodynamic potentials are related to the on-shell Euclidean 
action (times the temperature to cancel the periodicity of the Euclidean time direction). 
Given our previous discussion on the identification of the magnetization with the field J (see 
(3.5)), which leads to the following relation 

dM _ M 
~dl ~ T ' 

we can associate the Legendre transformation defining the Gibbs free energy with the Legendre 
transformation defining the action S in (2.12). By denoting the on-shell action 7 as / we have 8 



(3.13) 



-J 



-N?ir 2 T 4 



1 + e h 




1 + ^ + ^log 



rt+ri 



+ 0(ei 



(3.14) 

The regularization is performed by subtracting the T = background and the action is sup- 
plemented with a Gibbons-Hawking term. In [20, 21] the regularization of the free energy 
in the probe approximation was performed by the addition of counterterms. We consider 
the fact that we recover their results in the appropriate limit as a sign that the background 
subtraction method gives the correct answer. In particular, we do not need to worry about 
the presence of logarithmic divergences, cancelled by counterterms with explicit cutoff de- 
pendence, since these are temperature independent and the background subtraction cancels 
them completely. It is not difficult to check that, when subtracting the T = background, 
the contribution of that term vanishes up to order 1/r*. Now, it is not difficult to check that 
indeed 

= ~A ST ~ 2 BM ■ (3 - 15) 
Applying the standard thermodynamic relations 

'dg~ 

dT 



M 



(3.16) 



dG_\ 

b \9BJ T 

we confirm the previously obtained results in (3.2) and (3.6) respectively. The calculation 
of the Helmholtz free energy can be done in a similar fashion, but using the original action 
(2.10). This can be seen as the Legendre transformation of S e ff, which eliminates Hq and 
adds a term + BAA, after the proper renormalization. In this way we have 



T = Q + BM =U - sT => T = 
Once again we can check the thermodynamic relations 

B = 



-sT + ]-BM . 
4 2 



dF_\ 



dMJ, 



(3.17) 



(3.18) 



7 Notice that to obtain S e f f we have integrated by parts to get rid of second order differentials, introducing 
some boundary terms. These, in principle, are taken care of by the Gibbons-Hawking term and will not 
contribute to the final expression. We have checked that this is the case by calculating the on-shell action with 
and without these extra boundary terms. 

8 This calculation is detailed in appendix C. 
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where, to work at fixed magnetization, we have to specify how r m evolves with the tempera- 
ture. For that we look at equation (3.6), from where the following evolution follows 

irR 2 r h r m 



drr ri 



■1 + log 



,+r' 



(3.19) 



3.3 



Speed of sound 

Finally we analyze the speed of sound in the plasma with a magnetic field. Due to the 
anisotropy of the gravitational solution we will find that there are two normal directions in 
which the pressure waves propagate at different speeds. For a perturbation propagating in 
the direction of the magnetic field we have 

2 _dP\\ _ -(dG/dT) B _ s 

s <ll dU (dU/OT) B c V)B ' K ' 

where Cy,B is the heat capacity at fixed magnetic field. To calculate it we have to derive 
the internal energy with respect to the temperature, but we must take into account how the 
parameters and r m run with the energy scale. The case of the parameter is easy to 
understand from the profile for the dilaton and it follows that &T^h ~ ci, [15]. Since we work 
at first order in e^, the running of the coupling constant 9 - via the presence of factors of R 
in the definitions of the physical magnetic field (3.4) and the magnetization (3.6) - does not 



affect our results. At fixed magnetic field, since B 



we observe that 8t r m = 0, therefore 



Cv,B = 



dU 
df 



B 



3iV c 2 rI 
2vri? 6 



( 



+ 3 Vi 



\ 



1 + 



+ 0(e 



7 



(3.21) 



With this result at hand we can find readily the speed of sound in the direction of the magnetic 
field as 



B 



4 



1 



6 r\ 



+ 0(e 2 h ) 



1 + 



(3.22) 



which gives a lower speed of sound than the conformal result. For the speed of sound in the 
direction orthogonal to the magnetic field we obtain, using the chain rule 

2 dp ± — (dj-/dT) B {WhMmkMlB 

(dU/dT) B {BU/dT) B 



B 



s ^ dU 
which leads to 



Cv,B Cy 



B 



7e h r: 



6 ri 



+ 0{e 2 h ) 



1 + 



(3.23) 
(3.24) 



'Recall that e h ~ X h N f /N c 
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In particular we see that the presence of a magnetic field in our setup breaks conformal in- 



variance at first order in A^-^ 



even when the fundamental degrees of freedom we included 
are massless (in the absence of magnetic field the breaking of conformal invariance happens 
at order e\, see [15-17, 22]. This is one difference between the setup presented in this work 
and the results in the quenched approximation XhNj/N c — > 0, [21]. Although for thermody- 
namic quantities such as the entropy, the magnetization and the Gibbs free energy we obtain 
agreement with the results of that paper, in our setup the anisotropy sourced by the magnetic 
field is included, and this allows us to calculate the different speeds of sound, depending on 
the direction of the pressure wave, and obtain conformality-breaking results. 

For completeness we calculate here the heat capacity at constant magnetization, where 
we need to make use of equation (3.19) to work at fixed magnetization. The result is 



M 



) m 



3iVfr 



2vri? 6 



24 



l + ^C^ M +0{e 2 h ) 



(3.25) 



with 



/~icor 



+ 



1 + 



3r\ [rl + 



r h + r r 



log 



,+r 



Gr 2 h {ri + 3ri) 



Sri 



rf + r. 



) - 2 (3ri + ri) 



log 



rt + 



+ r 



(3.26) 



4 Stress-energy tensor with a magnetic field 

In this section we will calculate holographically the stress-energy (SE) tensor of the boundary 
field theory. As customary in the AdS/CFT context, we evaluate the Brown- York tensor at 
a cutoff r\ from the 5d action (2.14) 



2 5S 5d 



"7 hi- 



(K ij - /\' -'•' ) 



'A 



(4.1) 



'A 



where 7 is the induced metric at the r = r\ surface, where the indices i,j run, and K lJ 
is the extrinsic curvature. The Brown- York tensor diverges when the cutoff is taken to the 
boundary. To cancel this divergence we employ the same background subtraction as in the 
previous section, which allows us to read the temperature and magnetic field contribution to 
the SE tensor. 



4.1 Expectations 

Before presenting the actual calculation we will state what we expect the diagonal components 
of the SE tensor to be. The presence of the magnetic field sources an anisotropy in the medium, 
and therefore we will have a vev for the SE tensor of the field theory given by 



(Tj) = diag {-Sadm, P±,P±, P\\) , 



(4.2) 
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where Sadm is the enthalpy, as was shown in the previous section, P± the pressure in the 
directions transverse to the magnetic field, and jPIi the pressure along the direction of the 
magnetic field. When no magnetic field is present the two pressures coincide and are related 
to the Gibbs free energy 10 P± = P\\ = —Q. 

The question that is immediately risen is whether in our case P± = —Q or Pn = —Q 
-if any-, and if this is true what is the expression for the other pressure. Notice that the 
difference Ap = Pn — P± is a measure of the anisotropy of the medium, and therefore we 
expect it to be proportional to the magnetization (times the magnetic field) Ap ~ BM. 

The answer to this question is given by 

P\\=S , P± = -J r . (4.3) 

To understand why this is the case, we will follow a thermodynamic argument that can be 
found in a similar context in appendix C of [25]. In that paper the thermal N = 4 SYM 
plasma has an anisotropy sourced by a specific distribution of D7 branes along the horizon of 
the black brane, translated in a value for the axion \ = az > with a a constant. In the present 
case the D7 branes are extended along the radial direction of AdS, reaching the boundary and 
describing fundamental matter in the plasma, and the anisotropy appears by the presence of 
a magnetization of the fundamental. 

The key of the argument is to write the internal energy of the plasma as an extensive 
quantity 11 U = U{S,L x ,L y ,L z ,M), with S the extensive entropy, L XjV:Z the length of the 
sides of a box in which we have inserted our plasma and M the magnetization of the system. 
The energy and the entropy scale with the total volume of the box V3 = L x L y L z = f d?x, but 
the magnetization does not. This may seem strange at first sight, since one would expect the 
magnetic field to be an intrinsic quantity and the magnetization to be a density. One way to 
see the scaling is to realize that the magnetization is a vector in the z direction whereas the 
magnetic field is given by a 2-form B dx A dy. Therefore, to keep the magnetic field constant 
when we scale L x or L y , we should scale B accordingly. In the same way, the magnetization 
scales with L z . This suggests that it is more appropriate to talk about magnetic flux along the 
xy plane, B, and magnetization linear density along the z direction. These are the quantities 
that matter when considering a finite box is the presence of an external magnetic field. 

Therefore, comparing with the calculation in [25], all we need to do is to repeat the 
arguments in their appendix C with the identification a — > M and $ — > B -which we will 
not write explicitly here since it is nicely discussed in the referred paper-, and we are led to 
the result (4.3). From here, it is also straightforward to see that Ap = BM. 

The identities (4.3) can also be written as Gibbs-Duhem equations 

U + P\\ =sT + BM , U + P ± = sT. (4.4) 

10 In the absence of a magnetic field the Gibbs and Helmholtz free energies presented in (3.14) and (3.17) 
coincide, but the presence of non-trivial charge density and chemical potential would make a difference between 
the two. Actually, in a traditional nomenclature the Gibbs free energy should correspond to the thermodynamic 
potential at fixed chemical potential and zero magnetic field; we use the same name here by analogy. 
J Note that in the rest of the paper thermodynamic quantities are intensive! 
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4.2 Holographic calculation 

We proceed now to calculate the components of the vev of the SE tensor in the field theory. 
This is related to the Brown- York tensor (4.1) by 

{T 3 ) = v^rVe 9 , (4.5) 

where we have assumed that the expression (4.1) has been regularized before taking the 
r\ — y oo limit. Notice that strictly speaking this is a density since we are not integrating over 
the space. From the definition (4.1) we have 

(T* t ) = ^r 5 (l - !|) d r log (e-Wh-*W) , (4.6) 

(T x x ) = <Z\> = ^ (l - d r log (e-^h-^bbl) , (4.7) 

(T z z ) = ^r 5 (l - ^ a r log (e- 10 ^- 3 / 2 i» 2 62 ) . (4.8) 

Expressions (4.6)-(4.8) can be expanded in powers of using the solution described in ap- 
pendix B. With this we can write (r l j) = (r J j) + e/j(r l J ) 1 + 0{e\). 

At zeroth order in we have 6 = 1, since this function describes the anisotropy between 
the directions perpendicular to the magnetic field and the direction along the magnetic field, 
which is an order effect caused by the presence of fundamental matter. Therefore, at zeroth 
order in one obtains that (T x x ) = (T y y ) = (T z z ) and the Brown- York tensor is isotropic. 
Actually, at zeroth order the solution to the type IIB action is nothing but AdS§ x S 5 by 
construction, and we know already what the Brown- York tensor is going to be. The explicit 
calculation goes as 

(v /3 7^) odw = ^(3rl-^)diag(-3,l,l,l) , (4.9) 

K 5 

where the subindex div signs that the expression is divergent in the r\ — > oo limit and must 
be regularized. Once again, the regularization is achieved by background subtraction 




7 AU - ^1 - J r Jim M (v^ry^j , (4.10) 

where the factor in the square root matches the euclidean geometries at the cutoff. A straight- 
forward calculation gives 



( Tl j)o = ^4 diag (-3, 1, 1, 1) , (4.11) 

after use of k| = k\q/Vse- Considering now the observation made in (3.3), we can rewrite 
this expression as 

rri 

(r j ) = ^diag(-3,l,l,l) . (4.12) 
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This, of course, is just the AdS§ result, which gives an isotropic contribution. Notice that 
even when the entropy density and the temperature are sensitive to the magnetization of the 
plasma, their product cancels out factors coming from Ai to give the contribution to the SE 
tensor given above. 

We consider now the contribution due to the presence of fundamental matter at first 
order in e^, (T%)i. This term is given prior to regularization by 

(^=7^)^ = rA(r £; ^ [(6i - 3A; + T[) I 4x4 - Jading (0, 1, 1,0)] . (4.13) 

Once regularized we read the vev of the field theory SE tensor. In this case the (T t) compo- 
nent must coincide with the ADM mass calculation (see appendix D) given in the previous 
section. We have checked this explicitly by regularizing (4.13) and evaluating the expression 
one gets in terms of b±, Ai and Ti. This fact helps us to find the expressions for the SE 
tensor with the aim of the following two properties 

(V^'Ol,*, = (V=JT z z ) ldm , (4.14) 

(v^r^-)x^ = (V^TV y ) ldiv = (V=ir%) ldiv - . (4.15) 

Expression (4.14) tells us that the contribution at first order in for (T*t) an d {T z z ) 
coincide, and since we know that the time component is given by the ADM energy, which we 
already calculated, we get 

(T t t ) = -e ADM = ~8T+^BM {T z z )= l -sT+^BM = -G , (4.16) 

as announced. 

To evaluate the pressure in the transverse directions P± = {T x x ) we can make use of 
the relation (4.15). There are two equivalent ways to obtain the answer. The first and more 
obvious one is to evaluate the b± contribution in the r.h.s. of (4.15) and regularize. This can 
be seen to lead to 

{T\) = (T z z ) - BM = -T , (4.17) 

which is the expected result. Unfortunately, the evaluation makes use of the analytic -but 
somehow complicated- form of b\(r), and intermediate steps to arrive to this result imply 
writing down long, non-illuminating expressions. A second strategy would be to notice that 
the contributions to the regularized SE tensor from Ai and Ti vanish. This implies that 
the correction at order to (T x x ) is opposite in sign to the correction to {T z z ), giving once 
again (T x x ) = ^sT — ^BAi = —T. However, the explicit solution for Ai and Ti is more 
complicated that the one for b\, and intermediate expressions are again cumbersome equations 
which would lengthen this section without adding anything relevant. 

Of course, from the former arguments it follows that the anisotropic measure is given by 
Ap = BA4, in agreement with the thermodynamic argument of the previous subsection. 
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5 Energy loss in the magnetically anisotropic plasma 

In this section we will focus on calculating the energy loss of the partons as they propagate 
in an anisotropic plasma. 

An estimation of this influence is coming through the calculation of the jet quenching 
parameter q. In [26], using the eikonal approximation in the high energy limit, they presented 
a non-perturbative prescription for calculating q as the coefficient of I? in an almost light-like 
Wilson loop with dimensions L~ S> L. Following this prescription we will calculate q for our 
anisotropic backreacted background. 

Another estimate of the energy loss of a parton passing through a medium may come 
through the drag force calculation. This computation can be implemented in a holographic 
framework through a macroscopic string moving with constant velocity v. That string is 
attached to a probe flavor brane and dragged by a constant force / which keeps the velocity 
fixed. The drag coefficient //, measuring the energy loss, is calculated from the equation of 
motion 



where p is the parton momentum. Following the prescription of [5], we will calculate \x for 
our anisotropic backreacted background. 

5.1 Drag force 

In this section we perform the calculation of a second observable describing energy loss in 
the Quark-Gluon plasma by computing the drag force experienced by an infinitely massive 
quark propagating at constant velocity through an anisotropic plasma in constant magnetic 
field. In an anisotropic medium, the drag coefficient is not just a number but a matrix. This 
matrix is diagonal, fi = diag(/i x , fi y , fi z ) but with fi x = fi y ^ fi z , which means that the force 
and the momentum (or the velocity) of the quark will not be aligned in general. 

The external magnetic field plays a double role in this scenario. On one side it makes 
the plasma anisotropic, on the other side it stimulates synchrotron radiation of gluons which 
is an additional factor contributing to the energy loss of the moving quark. In our analysis 
we stabilize the classical trajectory of the quark by introducing an additional electric field 
perpendicular to the magnetic field, compensating the Lorentz force. In addition we add a 
drag electric force compensating the viscous force of the plasma. The only energy loss is due 
to the negative work exerted by the viscous force. 

We follow closely [27], where the isotropic analysis of [5] has been generalized to the case 
of anisotropic plasma. Another relevant papers are ref. [28], where heavy quark in external 
magnetic field has been studied, and ref. [29] where the study of the radiation of a quark in 
an anisotropic plasma is performed. 

On the gravity side the quark is described by a string propagating in the background 
(5.23) while the string action is given by 



/ = VP , 
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where g is the induced worldsheet metric and P[B] is the pullback of the Kalb-Rammond 
B-fidd. 

Physically, the electric forces needed to stabilize the trajectory are introduced by attach- 
ing one end of the string to a D7-brane and turning on a constant gauge filed Fmn = d^Am 
on the brane [27]. This results to the following boundary term 



S 



bdry 



drAfqd T X 



N 



Next we define 



n 



1 

2 
dC 



dTF MN X M d T X 



M 



d(d a x M ) ' 

From the variation of the boundary action one has 

n M |aE + (2ira')F M NdrX N = . 



(5.3) 
(5.4) 

(5.5) 



Equation (5.5) is the equation for the balance of the forces acting on the moving quark. 
One can expect that the contribution to Hm from the Nambu-Goto term of the action cor- 
responds to the viscous force of the plasma, while the contribution from the anti-symmetric 
part corresponds to the Lorentz force, this is confirmed by our analysis. 

We choose an ansatz in which the string does not move along the compact directions, 
while due to the rotational symmetry in the xy-plane, we can choose y = 0. However in order 
for this ansatz to be consistent we need to compensate the Lorentz force along the y direction. 
To clarify this we keep a general ansatz y = y(r, a) for a while. We fix the reparameterization 
invariance by identifying (t, r) = (r, a) and consider a string profile of the form 



[ut + x(r)] simp 



z = [ut + z(r)] cosy? , y = y(r,a) , 



(5.6) 



corresponding to a quark moving with velocity u in the xz-plane at an angle <p with the 
z-axis and with so far undetermined profile along y (eventually we will fix y = 0). Since the 
lagrangian does not depend explicitly on x , y and z we have: 
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n 
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_Gn 
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_ ^ 22 
y ~~~ 
+ u simp B 12 

1 



G t t x' + u 2 cos 2 ip G 33 (a/ - z') + G 2 2V (yx' - uy') sinp - B u y , (5.7) 



Gu y' + u sin 2 ip G\\ (uy' — x'y) + u cos 2 p> G 



33 (uy ~ z 



'y) 



c 



G 



33 



G t t z' - u 2 sin 2 ip Gu (x' - z') + G 2 2V {yz' 



uy 



cos ip , 



(5.8) 



(5.9) 



where ' denotes differentiation with respect to r. As one can see from equation (5.8), even 
if we set y = 0, Tl y has a non- vanishing contribution usimp B12. For the y component of 
equation (5.5) we obtain 

u sin ip H - (2ira')F 2 = , (5.10) 
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where we have used that Fmn has only electric components and that -B12I9S = Hq. It is clear 
that equation (5.10) represents the balance between the Lorentz force usinpHo acting on the 
quark and the electric field along the y component, Fq2, needed to cancel the Lorentz force. 
Therefore we fix Fq2 = u sin p Ho/(2ira') which enables us to set y = 0. The expressions for 
H x and Hy simplify to: 



n, 



G tt x' + u 2 cos 2 tp G33 (x' - z') 
r tp Gn (x' — z') 



Inverting (5.11) and (5.12) we have 

Gaa G rr u 2 N 2 



(x') 
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Gn G tt D - N X N Z 



CO 



f\2 



smip , 
COS if , 

Gn G rr U 2 N 2 Z 



G33 Gtt D - N X N Z 



with 



and 



n^n, 



(5.11) 
(5.12) 

(5.13) 

(5.14) 
(5.15) 

D = cot ip Gtt Di D 2 (5.16) 
u l Gn G33 sin 2ip and D2 = Gtt + u 2 (G33 cos 2 tp + Gn sin 2 ipj . 

N z = and D = N x N z = 



N x = H x Gtt cot ip + u 2 Gn cos (f (Tl z cos(p + H x simp) 
N z = U z G tt + u 2 G 33 cos <p (U z cos (p + U x sin <p) , 



I ,.2 
2 



The critical value of r is determined from the equations N 3 
and it is given by 



(l-u 2 ) 



2^/4 



4 2 2 
! , c r m « COS^ <y3 



h 



1-U 2 



ir 2 + log 2 (1 - u 2 ) + Li 2 



u 



1 



(5.17) 



For this critical value we have Tl x = uGnsvmp and Tl z = uG^ cos (p. With all these con- 
straints the denominator in (5.13) is always real and positive except at r c , where it vanishes. 
The numerators, at the critical point, also vanish and the functions x',z' are smooth and 
negative for < r < r s — > 00. The force that must be exerted on the quark in order to 
maintain its stationary motion (see [27] for a detailed explanation) is 



F 



2ira 



7 (n^iT) , 



(5.18) 



in terms of the quark's velocity u = u(sin tp, cos (p). For the specific case of anisotropic plasma 
under study we have, in the small magnetic field limit 



X h T 2 



l + \e h [2- log(l - u 2 ) + f {XjZ} ] 
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with 



/. = ^7jT 36-121og(l-n 2 ) - 6 f 1 + ^ 
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(5.20) 
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2 ) + Li 2 
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. (5.21) 



These expressions show the dependence of the magnitude of the force on the magnetic field, 
and also how the presence of the anisotropy affects the orientation of the field. This effect 
is not due to a Lorentz force, since it is compensated in equation (5.10). The directions in 
which the quark is moving and the force is pushing are related by 



where bi(r c ) is the correction to the function b in the ansatz for the metric, evaluated at the 
critical radius. 

Finally we would like to discuss the regime of validity of our analysis. In this section 
we introduce additional electric fields coupled to the charged quark, but we do not take into 
account the effect that they have on the SYM plasma. On the other hand we take into 
account the effect of the magnetic field on the plasma. This can be justified if the quark 
is moving sufficiently slow and the Lorentz force and the viscous force (both proportional 
to the velocity of the quark) can be compensated with electric fields weak relative to the 
magnetic field. Comparing to the discussion in section 2.4, for our solution to still be reliable 
we consider that the electric field enters, effectively, at order e 2 . This is the regime in which 
our analysis applies. 

5.2 Jet quenching 

In this section we will compute the jet quenching parameter q for our anisotropic plasma, 
extending the computation of [15]. 

We will follow the analysis of ref. [30], where jet quenching in anisotropic plasma have 
been studied. In our case the anisotropy is due to external magnetic field which triggers 
synchrotron radiation. Therefore one would expect the motion of the fundamental string 
to encode both the phenomena of jet quenching and the energy loss due to synchrotron 
radiation. This suggests that the prescription of refs. [26], [30] should be revised in the 
presence of external magnetic field. Technically this can be seen from the more general form 
of the fundamental string action (5.2). In particular the presence of the -B-field term affects 



arg F = (p + e h 




sin(2^) + 0{e\) 



(5.22) 
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the prescription for the holographic calculation of the wilson loop. More precisely the on- 
shell action evaluated on the world sheet corresponding to the wilson loop is no longer purely 
imaginary which is crucial part of the derivations in refs. [26], [30]. 

To circumvent this difficulty we will constraint ourselves to motion parallel to the mag- 
netic field. In this case the pull-back of the -B-field on the world sheet vanish and the motion 
of the string in not directly affected by the external magnetic field. Physically the Lorentz 
force acting on the quarks vanish. This suggests that in the holographic calculation the wilson 
loop is given by the area of the world sheet and technically there are no difficulties in applying 
the prescription of refs. [26], [30]. Note that the energy loss of the quark is still indirectly 
affected by the external magnetic field through the alternated properties of the SYM plasma. 
We will sketch the derivation referring to [30] for details 12 . 

We consider a parton moving along the direction parallel to the magnetic field, z, with 
the momentum broadening taking place in the transverse xy-plane. Due to the rotational 
symmetry in the transverse plane we can choose q to lay along the x-direction. 

In order to cover more general situations we consider the following class of metrics 

ds\ Q = G tt dt 2 + Gn dx 2 + G 22 dy 2 + G 33 dz 2 + G„ dr 2 + • • • (5.23) 
where the ellipses denote compact directions. After introducing light-cone coordinates 

z ± = (t ± z) , (5.24) 

we consider a rectangular Wilson loop with contour C. The expectation value of the Wilson 
loop is given by the extremum of the Nambu-Goto action for a string with endpoints tracing 
the contour C. We consider a quark moving along z~ and fix the reparameterization invariance 
by identifying (z~,r) = (t, a). Moreover we set z + = and specify the string embedding 
through the function x = x(r), subject to the boundary condition x(±£/2) = 0. The Nambu- 
Goto action for this configuration is then given by 



S = 2i ^ J dz^ dr \j\(G u + G 33 )(G rr + G u x' 2 ) , (5.25) 

where x' = dx/dr and the action is imaginary because the string worldsheet is spacelike. 
Since the action does not depend on x explicitly we obtain the equation of motion for x(r) 

ir 2 G 
= 4 Gn (G tt + G 33 ) - vr2 GIT ' (5 ' 26) 

The turning point for the string is at x' = and following the prescription of [26] and [33] 
we will work in the limit I — > 0, which corresponds to the limit tt x —> 0. Integrating (5.26) 
and taking the limit for ir x — > 0, we obtain the separation length between the endpoints of 
the string 

L = tt x 1 x + 0(tt 2 x ) with l x = ( * ^£ dr . (5.27) 

Jr h Gu yGtt + (-T33 



2 See [31, 32] as well. 
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For the computation of the jet quenching parameter we have to evaluate the action on shell 
and focus on the I? term after using (5.27). In way we have 



,5' 



ilr \f\ 11/ 



8V2 2ttR 2 1 x ' 

The prescription given in [26] and [33] for the jet quenching parameter is 

L~i 2 .1 2 



J2S 



cxp 



4y/2 



2ttR 2 l r 



(5.28) 



(5.29) 



We rewrite this expression in terms of field theory quantities, namely the 't Hooft coupling 
at the temperature scale, A^ = Ae*' 1 and the temperature, given by (3.1). We obtain a 
correction to the unflavored jet-quenching in the presence of a small magnetic field given by 



r(f) 



A,T 3 



1 + 77^ 



r 4 / 1 
2 + 7T + 3 - - M cor 
rf \ 6 



with 



vr (log 8 - vr) + 4^3 



5 7 

1 1 1 ---22-1 



+ O ( 

r h 



-1.99143 . 



+ 0(e 2 h ) 



(5.30) 



(5.31) 



Unfortunately we were not able to find a closed expression for the jet quenching parameter 
when generic magnetic field is present. 

Comparison to the non-magnetic case 

To understand whether the presence of the magnetic field enhances or reduces the energy loss 
parameterized by q with respect to a theory without magnetic field, we must compare the 
expression (5.30) obtained before between the two different theories. 
The result in (5.30) can be written as 



q0 [l + e h ^M + O(ei] 



where 



7r 3 / 2 r(!) r- , 



(5.32) 



(5.33) 



is the flavored result in the absence of a magnetic field [15], which receives a enhancement of 
the jet quenching with respect to the unflavored setup, both in a scheme where the number of 
degrees of freedom (entropy density) and temperature are kept fixed or in a scheme where the 
energy density and the force between external quarks are fixed, indicating the robustness of the 
correction. The presence of the magnetic field is given by the factor M. = (3/8 — M cor /48)/8 ~ 
5/12 when the quark is moving along the direction of the transverse field. 

To compare the magnetic and non-magnetic result we must state clearly under which 
conditions we are making this comparison. For example, if we choose to keep the entropy 
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density and the temperature of the field theories the same, therefore allowing us to compare 
the values of q per degree of freedom at fixed T, we observe that the parameters N c in the 
B = and B ^ cases are related (in the r m < approximation) by 



N, 



c,B 



N, 



c,B=0 



1 



4 r\ 



+ O (el) 



(5.34) 



1 /2 

This correction enters in the jet quenching parameter via y/\h ~ N c , giving 

„4 



QB 

4b=o 



1 + e h - 



h 



M 



+ 0(e 2 



h) 



(5.35) 



implying that the presence of a magnetic field enhances the jet quenching if the quarks are 
moving parallel to the magnetic field. 

We could have chosen to fix N c as well as the entropy density, allowing T to vary. In 
that case the temperatures in the presence and absence of a magnetic field are related by 



6 r\ 



QB=o 



l + e h ^ M-- +0(el 



(5.36) 



Therefore, in this scheme the presence of an anisotropy induced by the magnetic field reduces 
the jet quenching for motion parallel to magnetic field. 



6 Conclusions 



In this work we have presented a solution to the equations of motion of type IIB supergravity 
in the presence of a smeared set of Nf S> 1 D7 branes. This solution is perturbative in 
the backreaction parameter ~ ^hjf- The presence of a finite magnetic field sources an 
anisotropy in the solution which leaves a footprint in physical observables. 

We have studied the thermodynamics associated to the magnetically anisotropic solution 
at first order in backreaction. Results for thermodynamic quantities like the entropy, free 
energy or magnetization coincide with studies performed in the quenched approximation, 
in which the dynamics of the matter in the fundamental representation decouples from the 
dynamics of the adjoint degrees of freedom. However, the quenched setup fails to describe 
the anisotropy due to the magnetic field. Actually, this probe approximation is only valid 
at small magnetic fields (compared to the scale of the temperature), where the anisotropy is 
very mild and the component of the NSNS potential, B xy (r), is approximately constant (see 
figure 5). At larger magnetic fields the backreaction of the branes onto the geometry creates 
a non trivial profile of the B xy component, which is itself reflected in an anisotropy in the 
metric. 

We have presented expressions for the pressure of the plasma in the directions parallel and 
orthogonal to the magnetic field, that -not surprisingly- do not coincide for finite magnetic 
field. This has as a consequence that the speed of sound in the two normal directions of the 
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plasma do not coincide between them. Actually, from equations (3.22) and (3.24) we observe 
that both speeds of sound have a value lower than the conformal setup, signaling that scale 
invariance is broken by the magnetic field at first order in eh, even when we have massless 
D7 branes (in the absence of magnetic field this is not the case, even when a charge density 
is present in the setup). 

The breaking of the conformality means that we must not necessarily have a traceless 
stress-energy tensor. However, from direct inspection in section 4, we have (T^u) = at first 
order in backreaction. A word of caution is needed here. As discussed in some extent in the 
text, the 00 component of the stress-energy tensor corresponds to the magnetic enthalpy, and 
not the internal energy of the system. Therefore one must not conclude that the tracelessness 
of the s-e tensor implies an equation of state of the form IA = 2P± + P\\- In fact, given the 
relation between the magnetic enthalpy and the internal energy (3.10) we have 



The last consequence of the anisotropy we have studied is the implications of the magnetic 
field in the energy loss of a heavy quark moving through the plasma. The lack of explicit 
isotropy has as a consequence that the energy loss depends on the direction of movement of 
the quark. Furthermore, if the quark is charged it will feel a Lorentz force due to the presence 
of the magnetic field. 

An alternative approach to the jet quenching calculation we presented in this work is 
through fluctuation analysis, relating the jet quenching parameter to momentum broadening 
[34]. In this reference the jet quenching is related to the transport coefficient associated 
to Langevin diffusion inside the plasma. A similar analysis of the jet quenching parameter 
for a deformed Af = 4 SYM after introducing massless flavor branes in the Veneziano limit 
was presented in [35]. Since the background is analytic, but perturbative in the number 
of flavors, it is possible to obtain perturbative expressions for the jet quenching around the 
unquenched result. We believe that the same analysis could be very well extended in our case, 
with the only concern on the complexity of the solution after the inclusion of the magnetic 
field. Generically, it seems a straightforward computation that will circumvent the issue of 
generalizing the "standard" recipe of the jet quenching calculation. 

One interesting question that is raised immediately is what are the shear and bulk vis- 
cosities of the magnetic plasma. The shear viscosity is a tensorial quantity that has been 
seen in [36] not to satisfy the KSS value -q/s = 1 /An for the anisotropic plasma of [37] (this 
situation happens as well in condensed phases of holographic superfluids, as was proposed in 
[38] and checked explicitly in [39]). Here we have an anisotropy sourced by a 2- form instead 
of a scalar, and this may complicate the analysis of the perturbations to calculate the shear 
viscosity via a Kubo formula. Due to the lack of conformality at first order in e^ we expect 
that the bulk viscosity is non-zero as well, but proportional to BM. 



U = 2P ± + R, + BM . 
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A Equations of motion from the effective action 

The expression for the one-dimensional effective lagrangian is given by C e ff 

' Ql \ 2 77/ Ql 

S » ~ r - • T.2 l2 ei2 c6 . 1,2 -2 P 4 r-4 
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1 2 l e-*H*h ( e 2 *J> 2 b 2 \ 
Defining the following auxiliary (dimensionless) expressions 



we can write the equations of motion in the following compact way 
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It is straightforward to check that the above set of equations, together with (2.11), solve the 
full set of Einstein equations, provided the following "zero-energy" constraint is also satisfied 
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This constraint can be thought of as the aa component of the Einstein equations. Differenti- 
ating (A. 10) and using (2.11) and (A.4)-(A.9) we are getting zero, meaning that the system 
is not overdetermined. 
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B Analytic perturbative solution of the equations of motion 



The homogeneous solutions for the equations (2.24)-(2.25) are 
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(B.2) 

(B.3) 
(B.4) 
(B.5) 
(B.6) 



The solutions with integration constant K# \ are regular at the horizon whereas the ones with 
K\s,2 diverge logarithmically there. This situation is reversed at infinity, where the solutions 
with K\$, 2 tend to zero whereas the ones with Km can diverge or go to a constant. Also, in 
principle Qi/zftr /f^ — 1) has an imaginary part, but this is just P ly / 2 (2r 4 /r^ — 1) and can 
be absorbed in Ka,i- In the following, we will consider just the real part of <3i/2(2f 4 / r h — 1) 
and that constants ^a,i,2 are real. 

The particular solution for every equation can be found in the following way. Defining 
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a particular solution of the corresponding inhomogeneous differential equation is given by 
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with ^ ' the homogeneous solutions accompanied by the constants (12) and 



(B. 



(B. 



the Wronskian. With the symbol J r we denote an antiderivative, therefore no lower bound is 
considered (its addition amounts to a shift in the constants of integration for the homogeneous 
solutions). This method also works for Hi, but in that case 
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The solutions can be expressed in terms of the following expressions 
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which for the special cases n = 1, 5, 9 read 
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We express the solution^ in terms of the dimensionless parameters ct r , ot Th and 



a,. 



Solution for 01 

The solution for <j)\ which is regular at the horizon and vanishes at r = r* is 
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Whenever we calculate a physical quantity we will set the scale r* at the horizon = r^, 
indicating that we are working in the IR range of energies. 



Solution for b\ 

The solution for b\ which is regular at the horizon and vanishes at r = r s is 
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13 We do not include the expression for Ai since it is given in terms of integrals that cannot be evaluated 
analytically. Furthermore, we do not need its explicit form in the text. 
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The following limit of (B.15) will be useful in thermodynamic calculations 
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Solution for Ai 

The solution for Ai 14 which is regular at the horizon and vanishes at r = r s is 
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The following limit of (B.18) will be useful in thermodynamic calculations 
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*The value of the constant Ca ± is determined by requiring Ai(r s ) = 0. 
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Solution for Ti 



The solution for Tx is 
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The following limit of (B.21) will be useful in thermodynamic calculations 
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Solution for ifx 



The solution for Hi 16 is 
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The limiting cases — > and r m — )• can be obtained, recovering the results in [19] and [15] 
respectively. 



'The value of the constant Cr 1 is determined by requiring Ti(r s ) = 
3 The value of the constant Ch x is determined by requiring H\(r s ) = 



- 35 - 



C Calculation of the Gibbs free energy 



Applying the general recipe of [40] we identify the on-shell Euclidean action, divided by the 
inverse temperature, with the Gibbs free energy (in an ensemble where the magnetic field is 
kept fixed). The Euclidean action has contributions from both bulk and surface terms given 
by the following expressions 
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2k 10 



Cub da and J sur f 



V ^ 3 ru 



(C.l) 



where Cub is the wick rotated action (2.6) and in X sur f is the standard Gibbons-Hawking 
term. In order to generate £ e //> defined in (A.l), we need to integrate by parts and get rid 
of second order derivative terms. This in turn will lead us to consider boundary terms, on- 
shell, only of the form ip' a /i^ a , where ip a is a collective notation for the background functions 
(b, br, h, F, S, <&). While at r = r s those terms are canceled by the contribution of the 
Gibbons-Hawking term, they remain at the horizon. It turns out that the only non-vanishing 
term is coming from the function bx, with the following contribution 



d a b 



T 



In other words 



1- = ^bulk + Zsurf = P 7T-1T ( 4 r h 



C e s da 



Expanding in and using (2.20) we have 



(C2) 



(C3) 



Z = ?0 + e/i^DBl + e/i^bound + 0(e 2 h 
and more explicitly, changing the radial coordinate to r 
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1 



[5bi(r s ) + 5Ai(r fl ) - Ti(r s )] 



(C.4) 



(C.5) 



This action is infinite and should be regularized by subtracting the zero temperature on-shell 
Euclidean action. We take r s as the radial cut-off for the integrals, such that the finite and 
zero temperature geometries coincide. Since gu{ r s) ^ 9tt{ r s), we rescale the Euclidean time 
of the zero temperature solution in the following way [15] 



1/2 



Po = P 1 



(C.6) 
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where (3 is the period of the Euclidean time of the finite temperature solution. Doing so it is 
easy to prove that 



1 



P Nl tt 2 T 4 



1 + e h 
and consequently (3.14). 



1 / r 4 r 4 



+ 0(e 2 h ) 



(C.7) 



D ADM energy and Brown- York tensor 

In this appendix we show that the definition of the ADM energy used in (3.7) coincides with 
the calculation of the tt component of the boundary Brown- York tensor, given by Eby = 
— J d 3 x-y/— 'jki^jT 11 , where k is the unit-norm vector orthogonal to the t = constant surfaces, 
£ = dt is a Killing vector and 

2 6S5d (D.l) 



\/-7 Hij 

the Brown- York boundary tensor. The calculation is done by noticing that in the boundary 
metric, 7jj, &T( r s) enters only in the jtt component, then we can trade "ftt by bx in the 
calculation 17 



J} 



b T (r s ) 



On the other hand 
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Sb T (r s ) iino 

where we are using (A.l) to define £ e //. Using the equation of motion for 6y we can express 
the on-shell value of C e f f as a total derivative and perform the integral and the variation 



SS 5 d _ V 3 V S E 9L e ff 



Sb T (r s ) 

Therefore, from the BY definition 
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Now, the ADM mass is defined in (3.7) as 
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rv 1Q rv 1Q 

Therefore, as far as we implement the same regularization procedure, these two quantities 
coincide and, as we have argued in the main text, can be identified with the magnetic enthalpy 
in the field theory side. 



17 Notice that we express the radial dependence of fields in terms of the radial coordinates r or a indistinc- 
tively, depending on which one is more convenient for every step. 



-37- 



References 

[1] J. M. Maldacena, "The large N limit of superconformal field theories and supergravity," Adv. 
Theor. Math. Phys. 2, 231 (1998) [Int. J. Thcor. Phys. 38, 1113 (1999)] [arXiv:hcp-th/9711200]. 

[2] G. Policastro, D. T. Son and A. O. Starinets, "The Shear viscosity of strongly coupled N=4 
supcrsymmetric Yang-Mills plasma," Phys. Rev. Lett. 87, 081601 (2001) [hcp-th/0104066]. 

[3] P. M. Chesler and L. G. Yaffc, "Boost invariant flow, black hole formation, and 

far-from-cquilibrium dynamics in N = 4 supcrsymmetric Yang-Mills theory," Phys. Rev. D 82 
(2010) 026006 [arXiv:0906.4426 [hep-th]]. 

[4] V. Balasubramanian, A. Bernamonti, J. dc Boer, N. Copland, B. Craps, E. Keski-Vakkuri, 

B. Mullcr and A. Schafer et al, "Holographic Thermalization," Phys. Rev. D 84 (2011) 026010 
[arXiv:1103.2683 [hep-th]]. 

[5] C. P. Hcrzog, A. Karch, P. Kovtun, C. Kozcaz and L. G. Yaffc, "Energy loss of a heavy quark 
moving through N=4 supcrsymmetric Yang-Mills plasma," JHEP 0607, 013 (2006) 
[hep-th/0605158]. 

S. S. Gubser, "Drag force in AdS/CFT," Phys. Rev. D 74, 126005 (2006) [hcp-th/0605182]. 

[6] V. P. Gusynin, V. A. Miransky and I. A. Shovkovy, "Catalysis of dynamical flavor symmetry 
breaking by a magnetic field in (2+l)-dimcnsions," Phys. Rev. Lett. 73, 3499 (1994) 
[Erratum-ibid. 76, 1005 (1996)] [arXiv:hep-ph/9405262]. 

V. P. Gusynin, V. A. Miransky and I. A. Shovkovy, "Dimensional reduction and dynamical 
chiral symmetry breaking by a magnetic field in (3+l)-dimensions," Phys. Lett. B 349, 477 
(1995) [arXiv:hcp-ph/9412257]. 

D. K. Hong, Y. Kim and S. J. Sin, "RG analysis of magnetic catalysis in Dynamical symmetry 
breaking," Phys. Rev. D 54, 7879 (1996) [arXiv:hcp-th/9603157]. 

K. G. Klimenko, "Three-dimensional Gross-Neveu model in an external magnetic field," Theor. 
Math. Phys. 89, 1161 (1992) [Teor. Mat. Fiz. 89, 211 (1991)]. 

K. G. Klimenko, "Three-dimensional Gross-Neveu model at nonzero temperature and in an 
external magnetic field," Z. Phys. C 54, 323 (1992). 

K. G. Klimenko, "Three-dimensional Gross-Neveu model at nonzero temperature and in an 
external magnetic field," Thcor. Math. Phys. 90, 1 (1992) [Tcor. Mat. Fiz. 90, 3 (1992)]. 

[7] V. G. Filcv, C. V. Johnson, R. C. Rashkov and K. S. Viswanathan, "Flavoured large N gauge 
theory in an external magnetic field," JHEP 0710, 019 (2007) [arXiv:hcp-th/0701001]. 

[8] V. G. Filev and R. C. Raskov, "Magnetic catalysis of Chiral Symmetry Breaking. A Holographic 
Prospective," Adv. High Energy Phys. 2010, 473206 (2010) [arXiv: 1010. 0444 [hep-th]]. 

[9] V. G. Filev, "Criticality, scaling and chiral symmetry breaking in external magnetic field," 
JHEP 0804, 088 (2008) [arXiv:0706.3811 [hep-th]]. 

J. Erdmenger, R. Meyer and J. P. Shock, "AdS/CFT with flavour in electric and magnetic 
Kalb-Ramond fields," JHEP 0712, 091 (2007) [arXiv:0709.1551 [hep-th]]. 
J. Erdmenger, N. Evans, I. Kirsch and E. Threlfall, "Mesons in Gauge/Gravity Duals - A 
Review," Eur. Phys. J. A 35, 81 (2008) [arXiv:071 1.4467 [hep-th]]. 

C. V. Johnson and A. Kundu, "External Fields and Chiral Symmetry Breaking in the 
Sakai-Sugimoto Model," JHEP 0812, 053 (2008) [arXiv:0803.0038 [hep-th]]. 



- 38 - 



A. V. Zayakin, "QCD Vacuum Properties in a Magnetic Field from AdS/CFT: Chiral 

Condensate and Goldstonc Mass," JHEP 0807 (2008) 116 [arXiv:0807.2917 [hep-th]]. 

V. G. Filev, C. V. Johnson and J. P. Shock, "Universal Holographic Chiral Dynamics in an 

External Magnetic Field," JHEP 0908, 013 (2009) [arXiv:0903.5345 [hep-th]]. 

V. G. Filev, "Hot Defect Superconformal Field Theory in an External Magnetic Field," JHEP 

0911, 123 (2009) [arXiv:0910.0554 [hep-th]]. 

N. Evans, A. Gcbaucr, K. Y. Kim and M. Magou, "Holographic Description of the Phase 
Diagram of a Chiral Symmetry Breaking Gauge Theory," JHEP 1003, 132 (2010) 
[arXiv:1002.1885 [hep-th]]. 

K. Jensen, A. Karch and E. G. Thompson, "A Holographic Quantum Critical Point at Finite 
Magnetic Field and Finite Density," JHEP 1005, 015 (2010) [arXiv:1002.2447 [hep-th]]. 
K. Jensen, A. Karch, D. T. Son and E. G. Thompson, "Holographic 
Berezinskii-Kosterlitz-Thouless Transitions," Phys. Rev. Lett. 105, 041601 (2010) 
[arXiv:1002.3159 [hep-th]]. 

N. Evans, A. Gebauer, K. Y. Kim and M. Magou, "Phase diagram of the D3/D5 system in a 
magnetic field and a BKT transition," Phys. Lett. B 698, 91 (2011) [arXiv: 1003.2694 [hep-th]]. 
N. Evans, T. Kalaydzhyan, K. -y. Kim and I. Kirsch, "Non-equilibrium physics at a holographic 
chiral phase transition," JHEP 1101 (2011) 050 [arXiv:101 1.2519 [hep-th]]. 
N. Evans, A. Gebauer and K. Y. Kim, "E, B, fi, T Phase Structure of the D3/D7 Holographic 
Dual," JHEP 1105, 067 (2011) [arXiv:1103.5627 [hep-th]]. 

S. Bolognesi and D. Tong, "Magnetic catalysis in AdS4," arXiv:1110.5902 [hep-th]. 
M. S. Alam, V. S. Kaplunovsky and A. Kundu, "Chiral Symmetry Breaking and External Fields 
in the Kuperstein-Sonnenschein Model," JHEP 1204, 111 (2012) [arXiv:1202.3488 [hep-th]]. 
S. Bolognesi, J. N. Laia, D. Tong and K. Wong, "A Gaplcss Hard Wall: magnetic catalysis in 
Bulk and Boundary," arXiv: 1204.6029 [hep-th]. 

[10] A. Karch and E. Katz, "Adding flavor to AdS/CFT," JHEP 0206, 043 (2002) 
[arXiv:hep-th/0205236]. 

[11] F. Bigazzi, R. Casero, A. L. Cotrone, E. Kiritsis and A. Paredes, "Non-critical holography and 
four-dimensional CFT's with fundamentals," JHEP 0510 (2005) 012 [hep-th/0505140]. 

[12] C. Nunez, A. Paredes and A. V. Ramallo, "Unquenched flavor in the gauge/gravity 

correspondence," Adv. High Energy Phys. 2010, 196714 (2010) [arXiv: 1002. 1088 [hep-th]]. 

[13] D. Arean, E. Conde, A. V. Ramallo and D. Zoakos, "Holographic duals of SQCD models in low 
dimensions," JHEP 1006, 095 (2010). [arXiv:1004.4212 [hep-th]]. 

J. Gaillard, D. Martelli, C. Nunez and I. Papadimitriou, "The warped, resolved, deformed 
conifold gets flavoured," Nucl. Phys. B 843, 1 (2011) [arXiv: 1004.4638 [hep-th]]. 
J. Schmude, "Comments on the distinction between color- and flavor-branes and new D3-D7 
solutions with eight supercharges," [arXiv:1007.1201 [hep-th]]. 

E. Conde, J. Gaillard, "Kutasov-like duality from D5-brancs wrapping hyperbolic cycles," Nucl. 
Phys. B848, 431-473 (2011). [arXiv:1011.1451 [hep-th]]. 

E. Caceres, C. Nunez, L. A. Pando-Zayas, "Heating up the Baryonic Branch with U-duality: A 
Unified picture of conifold black holes," JHEP 1103, 054 (2011). [arXiv:1101.4123 [hep-th]]. 

D. Elander, J. Gaillard, C. Nunez, M. Piai, "Towards multi-scale dynamics on the baryonic 
branch of Klcbanov-Strassler," JHEP 1107, 056 (2011) [arXiv:1104.3963 [hep-th]]. 

E. Conde and A. V. Ramallo, "On the gravity dual of Chern-Simons-matter theories with 



- 39 - 



unquenched flavor," JHEP 1107, 099 (2011) [arXiv:1105.6045 [hep-th]]. 

E. Condc, J. Gaillard and A. V. Ramallo, "On the holographic dual of N=l SQCD with 
massive flavors," JHEP 1110, 023 (2011) [arXiv: 1107.3803 [hep-th]]. 

A. Barranco, E. Pallante and J. G. Russo, "N=l SQCD-likc theories with Nf massive flavors 
from AdS/CFT and beta functions," JHEP 1109, 086 (2011) [axXiv:l 107.4002 [hep-th]]. 

[14] F. Benini, F. Canoura, S. Cremoncsi, C. Nunez and A. V. Ramallo, "Unquenched flavors in the 
Klebanov-Witten model," JHEP 0702, 090 (2007) [hep-th/0612118]. 

[15] F. Bigazzi, A. L. Cotrone, J. Mas, A. Paredcs, A. V. Ramallo and J. Tarrio, "D3-D7 
Quark-Gluon Plasmas," JHEP 0911, 117 (2009) [arXiv:0909.2865 [hep-th]]. 

[16] F. Bigazzi, A. L. Cotrone, J. Tarrio, "Hydrodynamics of fundamental matter," JHEP 1002, 083 
(2010) [arXiv:0912.3256 [hep-th]]. 

F. Bigazzi, A. L. Cotrone, "An elementary stringy estimate of transport coefficients of large 
temperature QCD," JHEP 1008, 128 (2010) [arXiv: 1006.4634 [hep-ph]]. 

[17] F. Bigazzi, A. L. Cotrone, J. Mas, D. Mayerson and J. Tarrio, "D3-D7 Quark-Gluon Plasmas at 
Finite Baryon Density," JHEP 1104, 060 (2011) [arXiv:1101.3560 [hep-th]]. 

[18] F. Bigazzi, A. L. Cotrone, J. Mas, D. Mayerson and J. Tarrio, "Holographic Duals of Quark 
Gluon Plasmas with Unquenched Flavors," arXiv:1110.1744 [hep-th]. 

[19] V. G. Filev and D. Zoakos, "Towards Unquenched Holographic magnetic catalysis," JHEP 
1108, 022 (2011) [arXiv:1106.1330 [hep-th]]. 

[20] J. Erdmcngcr, V. G. Filev and D. Zoakos, "Magnetic catalysis with Massive Dynamical 
Flavours," arXiv:1112.4807 [hep-th]. 

[21] T. Albash, V. G. Filev, C. V. Johnson and A. Kundu, 'Finite temperature large N gauge theory 
with quarks in an external magnetic field," JHEP 0807, 080 (2008) [arXiv:0709.1547 [hep-th]]. 

[22] D. Mateos, R. C. Myers and R. M. Thomson, "Thermodynamics of the brane," JHEP 0705, 
067 (2007) [hep-th/0701132]. 

[23] F. Bigazzi, A. L. Cotrone and A. Paredes, "Klebanov-Witten theory with massive dynamical 
flavors," JHEP 0809, 048 (2008) [arXiv:0807.0298 [hep-th]]. 

[24] A. L. Cotrone and J. Tarrio, "Consistent reduction of charged D3-D7 systems," 
arXiv: 1207.6703 [hep-th]. 

[25] D. Mateos and D. Trancanelli, "Thermodynamics and Instabilities of a Strongly Coupled 
Anisotropic Plasma," JHEP 1107 (2011) 054 [arXiv:1106.1637 [hep-th]]. 

[26] H. Liu, K. Rajagopal and U. A. Wiedemann, "Calculating the jet quenching parameter from 
AdS/CFT," Phys. Rev. Lett. 97, 182301 (2006) [hcp-ph/0605178]. 

[27] M. Chcrnicoff, D. Fernandez, D. Mateos and D. Trancanelli, "Drag force in a strongly coupled 
anisotropic plasma," arXiv:1202.3696 [hep-th]. 

[28] E. Kiritsis and G. Pavlopoulos, "Heavy quarks in a magnetic field," JHEP 1204, 096 (2012) 
[arXiv:1111.0314 [hep-th]]. 

[29] K. B. Fadafan and H. Soltanpanahi, "Energy loss in a strongly coupled anisotropic plasma," 
arXiv: 1206.2271 [hep-th]. 



- 40 - 



[30] M. Chernicoff, D. Fernandez, D. Mateos and D. Trancanelli, "Jet quenching in a strongly 
coupled anisotropic plasma," arXiv:1203.0561 [hcp-th]. 

[31] R. Baier and Y. Mehtar-Tani, "Jet quenching and broadening: The Transport coefficient q-hat 
in an anisotropic plasma," Phys. Rev. C 78, 064906 (2008) [arXiv:0806.0954 [hep-ph]]. 

[32] D. Giataganas, "Probing strongly coupled anisotropic plasma," arXiv:1202.4436 [hcp-th]. 

[33] H. Liu, K. Rajagopal and U. A. Wiedemann, "Wilson loops in heavy ion collisions and their 
calculation in AdS/CFT," JHEP 0703, 066 (2007) [hcp-ph/0612168]. 

[34] S. S. Gubser, "Momentum fluctuations of heavy quarks in the gauge-string duality," Nucl. 
Phys. B 790 (2008) 175 [hep-th/0612143]. 

[35] A. Magana, J. Mas, L. Mazzanti and J. Tarrio, "Probes on D3-D7 Quark-Gluon Plasmas," 
JHEP 1207 (2012) 058 [arXiv:1205.6176 [hcp-th]]. 

[36] A. Rebhan and D. Steincdcr, "Violation of the Holographic Viscosity Bound in a Strongly 
Coupled Anisotropic Plasma," Phys. Rev. Lett. 108 (2012) 021601 [arXiv:1110.6825 [hep-th]]. 

[37] D. Mateos, D. Trancanelli, "The anisotropic N=4 super Yang-Mills plasma and its 
instabilities," Phys. Rev. Lett. 107 (2011) 101601 [arXiv:1105.3472 [hep-th]]. 

[38] M. Natsuume and M. Ohta, "The Shear viscosity of holographic superfluids," Prog. Theor. 
Phys. 124 (2010) 931 [arXiv: 1008.4142 [hcp-th]]. 

[39] J. Erdmenger, P. Kerner, H. Zeller, "Non-universal shear viscosity from Einstein gravity," Phys. 
Lett. B699 (2011) 301-304 [arXiv:1011.5912 [hcp-th]]. 

[40] S. W. Hawking and G. T. Horowitz, "The Gravitational Hamiltonian, action, entropy and 
surface terms," Class. Quant. Grav. 13, 1487 (1996) [arXiv:gr-qc/9501014]. 



- 41 - 



